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Abstract. The strong lensing events that are observed in compact clusters of galaxies can, both statistically and individually, 
return important clues about the structural properties of the most massive structures in the Universe. Substantial work is ongoing 
in order to understand the degree of similarity between the lensing cluster population and the population of clusters as a whole, 
with members of the former being likely more massive, compact, and substructured than members of the latter In this work 
we exploit synthetic clusters extracted from the MareNostrum Universe cosmological simulation in order to estimate the 
correlation between the strong lensing efficiency and other bulk properties of lensing clusters, such as the virial mass and the 
bolometric X-ray luminosity. We found that a positive correlation exist between all these quantities, with the substantial scatter 
being smaller for the luminosity-cross section relation. We additionally used the relation between the lensing efficiency and the 
virial mass in order to construct a synthetic optical depth that agrees well with the true one, while being extremely faster to be 
evaluated. We finally estimated what fraction of the total giant arc abundance is recovered when galaxy clusters are selected 
according to their dynamical activity or their X-ray luminosity. Our results show that there is a high probability for high-redshift 
strong lensing clusters to be substantially far away from dynamical equilibrium, and that 30 - 40% of the total amount of giant 
arcs are lost if looking only at very X-ray luminous objects. 



1. Introduction 

Explaining the process of cosmic structure formation is one 
of the major successes of the standard model for cosmology. 
Accordingly, tiny dark matter density fluctuations produced 
during the inflationary era grew up due to gravitational insta- 
bility, eventually detaching from the overall expansion of the 
Universe and collapsing into bound structures. The observed 
statistical properties of these objects can help to understand the 

I large scale matter distribution, and the way in which this has 
evolved during the history of the Universe. 

Gravitational lensing is one of the most powerful tools for 

' studying the formation of cosmic structures in general, and the 
mass assembly history of the most massive ones, galaxy clus- 
ters, in particular. Weak gravitational lensing can be used to 
perform non -parametric reconstruction of the outer dark matter 
halo profile (Kaise r et alJl995l : lBartel mann & Schneider 2001; 
LOguri et al. 2009) , while strong lensing allow to put constraints 
on the inner slope (iRichard et al.ll2010r) . that is one of the main 
theoretical predictions of the cold dark matter paradigm. The 
statistics of strong lensing images also is a particularly interest- 
ing tool for the study of structure formation, since it is sensitive 



to both the abundance of galaxy clusters and their average in- 
ner structure. It is hence expected to give important hints on the 
assembly history o f cosmic structures as w ell as on the back- 
ground cosmology JBartelmann et alj|2003b . 



Much work has been devoted in the past decade to the grav- 
itational arc statistics, that is the abundance of strongly elon- 
gated images that are produced by galaxy clusters when act- 
ing as gravitational lenses on background galaxies. Early re- 
sults pointed toward a deficit of observed arcs with respect to 
the theoretical predic t ions for a standar d ACDM co smology 



(iLe Fevre et al.l Il994c iBartelmann et al 
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dFedeli et al.l 120061 : iFedeli & Barteknannl l2007h seemed able 
to cur e. Whether this discrepancy survives nowadays is un- 
clear (IFedeli et alJ 12008) . however this issue propelled sub- 
stantial effort in order to better understand the cluster popu- 
lation responsible for strong lensing events, and to what de- 
gree this ensenible d i ffers from the overall c luster population 
(IHennawi et al.ll2007l : iMeneghetti et al.ll2010l) . Moreover, con- 
siderable work has been done on the observational side as well, 
in order to explore the high-redshift tail of the cluster lens pop- 
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ulation JGladders et alj|2003l; IZaritskv & Gonzalej|2003l) . and 
to pave the way toward machine-based gravitatio nal arc finding 
dLenzen et al]l2004l:[Seidel & Bartelmannll2007h . 

The dependence of the cluster lensing efficiency on the 
inner structure, dynamical state and merger activity of indi- 
vidual objects has bee n thoroughly studied ( iTorri etalj 12004 : 
Meneghetti et al.l2007h . An important issue that has been how- 
ever scarcely investigated about strong lensing cluster is how 
the lensing efficiency correlates with bulk cluster properties, 
such as mass or X-ray luminosity. This aspect is extremely in- 
teresting, since it allows to better quantify the statistical prop- 
erties of the cluster lens population, and in principle would 
allow to directly map the abundance of galaxy clusters into 
the arc abundance on the sky without intermediate calcula- 
tions. In this work we tackle this issue making use of the large 
MareNostrum Universe cosmological simulation. This sim- 
ulation contains adiabatic gas in addition to dark matter, and 
hence allows to characterize the cluster lens population taking 
into account the role of baryonic physics on structure forma- 
tion. We additionally estimated the optical depth of the whole 
cluster population, which allows to compute the total number 
of arcs expected in the sky, and how this depends on selection 
biases such as the dynamical state or the X-ray luminosity. 

The rest of this paper is organized as follows. In Section|2] 
we summarize the properties of the MareNostrum Universe 
simulation, and the calculations relevant for strong lensing 
of simulated clusters. F urther details on this can be found in 
Meneghetti et al.l ( 120101) . In section [3] we summarize the red- 
shift distributions of the strong lensing cluster population. In 
Section |4] we report our results and in Section |5] we summa- 
rize our conclusions. In the remainder of this work, we refer to 
Meneghetti et all (l2010l) as to Paper I. 



2. Lensing in the MareNostrum Universe simulation 
The MareNostrum Universe ( Gottl ober & YepesI 120071) is a 



large scale non-radiative SPH cosmological simula tion per- 
forme d with the publicly available code Gadget2 dSpringel 
2005b . The simulation box has a comoving side of 500 Mpc 
/!"' and contains 2 x 1024^ particles, equally subdivided in 
number between dark matter and baryons. The mass of each 
dark matter particle equals 8.24 x IO^Mq/z"', and that of each 
gas particle, for which only adiabatic physics is implemented, 
is 1.45 X IO^Mq/;"'. The cosmological parameters of the sim- 
ulation are in agreement with the WMAP-1 year data release 
(iSpergel et al.i 2003) . namely Qni,o = 0.3, Qa_o = 0.7 and 
CTg = 0.9 with a scale invariant primordial power spectrum. 
The baryon density parameter is set to Qb,o - 0.04 5. How this 
parameter set deal with the later W MAP releases ( Spergel et al.l 
20071: iKomatsu et alj|2009l 1201 Ol) is discussed in Section |5j In 
the following we describe only the main issues of our strong 
lensing analysis, deferring further detail to Paper I. 

Bound structures within the simulation at each redshift 
snaps hot were identified with a FOF algorithm ( K lvpin et al 



19991) with a basic linking length of 0.17 times the mean in- 
terparticle distance. All the FOF groups with mass larger than 
IO'-^Mq/i"' were then stored into sub-boxes of cubic shape with 
side length 5 Mpc /i ' for the subsequent lensing analysis. For 



each structure, the three-dimensional density field was com- 
puted on a regular grid inside the relative sub-box, and then 
projected along the three orthogonal coordinate axes. Bundles 
of lig ht rays were than traced through the central part of each 
plane (iMeneghetti et al.ll2007h . evaluating for each of them the 
deflection angle as the sum of the contributions given by each 
surface density cell. The final deflection angle map has a side 
length of 1 .5 Mpc /i ' , which allows to fairly capture the details 
of the critical line structure (if present). 

A preliminary analysis, performed with low resolution de- 
flection angle maps has shown that 49, 366 clusters were able 
of producing critical lines for sources at zs =2 in at least 
one projection. For each of the projections of these clusters we 
computed high resolution deflection angle maps for the subse- 
quent strong lensing analysis. We evaluated the cross sections 
(Tk^^ for arcs with length-to - width ratio^ > Rp by adopting the 
semi-analytic algorithm of lFedeli et alJ (l2006r) . Such algorithm 
reproduces well the results of fully numerical lensing simula- 
tions, while being substantially faster. Moreover, being based 
only on the properties of the deflection angle map, it is equally 
applicable to analytic (F edeli & Bartelmannll2007l) and to nu- 
merical lenses. For such computations the equivalent angular 
size of sources was kept fixed to 0.5 arcsec, and the ellipticity 
thereof was randomly drawn in the interval [0, 0.5]. Given the 
extremely large number of projections that we considered, we 
were forced to examine only one single length-to-width thresh- 
old, Rq = 7.5, and one individual source redshift, that we set 
to zs = 2. These choices are rather common in strong lens- 
ing statistics studies, with Rq - 7.5 being a fair compromise 
between having good number statistics and strongly distorted 
images, and zs = 2 being typical for many real strong lensin g 
clusters (see, e.g.. lEliasd6ttir et al.l2007l : lLimousin et al.l2008l) . 
Of the 148, 098 (= 3 x 49, 366) projections that we have an- 
alyzed, only 11,347 produced a non-vanishing cross section. 
This number is however by far the largest ever considered for 
theoretical strong lensing studies to date. 

The total number of arcs in the sky with R > Rq that are 
produced by galaxy clusters that lens sources at a given red- 
shift Zs, and that are visible below some limiting magnitude 
cut OTiim, is given by A6?o(mii,n) = 4;rn(miim,zs)T^o(2s), where 
n('«iim, Zs) is the number of sources that, in the unit solid angle, 
have redshift zs and magnitude below the limit miim. Note that, 
in order for this estimate to be self-consistent, the number den- 
sity has to be always corrected for t he lensing magnifi cation 
bias ('Bartelma nn & SchneideJl200ll: iFedeliet alJl200 8V This 
bias has the effect of increasing the angular density of visible 
sources if the intrinsic number counts are steep enough, while 
decreasing it if the number counts are relatively shallow. 

The optical depth for long and thin arcs produced by 
sources at redshift zs is instead defined as 
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N(M, zl) , ° ,, , dMdzL, 



(1) 



where N(M, z) represents the total number of structures present 
in the unit redshift around z and in the unit mass around M. 
The function D\(z) is the angular diameter distance out to red- 
shift z, while (Tr„ is the average cross section of each cluster 
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Fig. 1. The redshift distributions of strong lensing clusters in 
the MareNostrum Universe. The black histogram considers all 
clusters displaying critical lines, the red one those with a non- 
vanishing cross section for giant arcs, and the blue one those 
expected to produce at least one giant arc on average. The thick 
solid green curve represents the lensing efficiency, given by the 
normalized lensing distance, while the dashed one represents 
the lensing efficiency weighted for the cluster abundance (see 
the text for details). 

in the MareNostrum Universe. This average is taken over the 
three orthogonal projections (see above). Please note that this 
is equivalent to the correct procedure of computing three differ- 
ent optical depths, using a different cluster projection for each 
one, and then average the result. 

When only a finite discrete set of masses is available, as is 
the case here and in all the numerical or semi-analytic studies, 
then the Eq. ([1]) reduces to 



tr„(zs) 



-m 



RoJ 



,rf4;rDi(zs) 



NiM,ZL)dM 
Jm, 



dzL, (2) 



where £ is the number of strong lensing clusters at each redshift 
snapshot, the masses have to be sorted from the smallest to the 
largest at each redshift step, M, < M,+i, and the quantity o"^ ^ 
is defined as 



O"fio,'- = 9 [o'RoW'Z) + CrR„(M,+ i,z)] . 



(3) 



In practice, this approach is equivalent to assigning the aver- 
age cross section of clusters with masses M,- and M,+i to all 
structures with mass between M, and M,+i. The optical depth 
presented in Eq. ([TJ will be a central quantity in the rest of our 
study and, as mentioned above, we shall set Rq - 7.5 hence- 
forth. 



3. The population of strong lensing clusters 

We begin by reporting the distribution in redshift of the galaxy 
clusters that are capable of strongly lens background sources 
in the MareNostrum Universe. As done in Paper I we selected 
three subsamples of such clusters: /) clusters that are capable of 
producing critical curves, //) clusters that have a non-vanishing 
cross section, and ///) clusters that, given suitable background 
source number counts, are expected to produce at least one 
giant gravitational arc on average. The cross section of this 
last class of objects should be around 0-75 ~ 10"^ Mpc^/z"^. 
Obviously, the statements above refer to capabilities of at least 
one projection of the selected object. The strength of gravita- 
tional lenses increases when moving from type /) to type Hi) 
since, e.g., the presence of a caustic does not imply a non- 
vanishing cross section, and the bulk of cross section values 
is actually lower than 10"^ Mpc^/i"^ (see Paper I). The redshift 
distributions for these three categories of lensing clusters are 
shown in Figure [1] 

The thick green curve in Figure [T] represents the lensing 
efficiency s{zl), that we simply defined as the effective lensing 
distance normalized such as to have unit integral between zl - 
and Zl = zs- In other words, the lensing efficiency is defined 



fi(ZL) = 



1 Da(zl)Da(zl,zs) 



T 



Da(zs) 



(4) 



where Da(zuZ2) represents the angular diameter distance be- 
tween redshifts zi and zi, and Da{z) = £>a(0, z), while 



T = 



1 



Da(zs) 



Jo 



A(Zh)DAiZh,Zs)dZh- 



(5) 



Since the effective lensing distance defines the criti- 
cal surface density, which is responsible for the for- 
mation of critical curves ( Subramanian & Cowlind Il986t 



Padmanabhan & SubramanianI Il988h . it is expected that the 
redshift trend of e{zL) at least resembles the distribution of 
strong lensing clusters. 

This is actually the case, with b(zl) dropping to zero both 
at Zl - and at zl - zs, as do the simulated distributions. 
However the redshift distributions of strong lensing clusters 
are substantially more peaked than the lensing efficiency de- 
fined above, and the latter displays a substantial high-redshift 
tail that is absent in the former. This fact can be understood re- 
calling that the lensing efficiency does not contain any informa- 
tion about the evolution of the cluster mass function. Massive 
clusters are progressively rarer with increasing redshift, im- 
plying that the redshift distribution of strong lensing clusters 
should also experience an extra drop as compared to that due 
to pure geometric suppression. This can be verified by addi- 
tionally weighting the lensing efficiency defined above with the 
number density of objects that at a given redshift zl have mass 
larger than some fixed threshold value. This is shown by the 
dashed green curve in FigurelT] where we adopted IO'^Mq/z"' 
as the lower mass limit. This roughly corresponds to the min- 
imum mass of strong lensing clusters in the redshift range 
0.2 < Zl ^ 1 ■ With this extra weighting the high-redshift tail of 
the simulated distribution of strong lensing clusters is correctly 
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Fig. 2. The correlation between the logarithm of the virial mass (x - log(M h /Mq)) and the logarithm of the average cross section 
for giant arcs (y = log(o"7.5 /i^/Mpc')) at several different redshift snapshots of the MareNostrum Universe. The red lines show 
the best linear fit, while the dashed green lines are the expected trend for a SIS lens model (see text for details). 



captured, as expected. There remain a disagreement at very low 
redshifts, due to the fact that very close to the observer the min- 
imum mass capable of producing strong lensing grows rapidly 
well above 10 '^Mo/i-'. 



Next, let us examine the differences between the redshift 
distributions of the three categories of strong lensing clusters 
that we have defined. Going from subsample to Hi), that is 
increasing the strength of lenses, we see that the distribution 
shrinks, with the high-redshift tail being progressively removed 
and the height of the intermediate-redshift bins increasing. This 
is an obvious consequence of the fact that the more efficient 
cluster lenses are expected to be found far away both from the 
observer and the source plane. As a consequence, no clusters 
belonging to sample ///) are found beyond zl ~ 0.9, while clus- 
ters merely producing critical lines extend up to zl ~ 1 -4. 



4. Results 

4.1. Lensing efficiency scaling relations 

4.1.1. Mass 

As a next point we studied the correlation between the virial 
mass of galaxy clusters in the MareNostrum Universe and 
their strong lensing efficiency, represented by the cross section 
for giant arcs. As done before, since for each model cluster we 
computed the cross section along three different (orthogonal) 
projections, we quantify the lensing efficiency of each individ- 
ual cluster as the average 0-75 over the three projections. In 
order to simplify the notation, we introduce two new variables, 
X = logiMh/M^) andy = log(a-7.5 h^/Mpc^). 

In Figure |2] we show the correlation between the variables 
X and y for the clusters of the MareNostrum Universe that are 
contained within various redshift snapshots. For each panel we 
also show the best linear fit, in the form y - a(zL)x + b(zL)- 
We also tried a higher order fit to the points in Figure |2l but the 
result is not significantly different. We explicitly included the 
lens redshift dependence of the two best fit parameters a and b, 
since in general their value will differ for different simulation 
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Fig. 3. The distribution of the logarithm of cross sections for giant arcs (y - log(o"7.5 h^/Mpc^)) around the best fit at fixed 
logarithmic mass (x = log{M h/Mo)). The green dashed curves show the best Gaussian fit, while the red solid curves show the 
best first-order correction to the Gaussian fit (see the text for details). 



snapshots. The redshift evolution of the best fit parameters is 
discussed in Section \42\ As expected, an obvious correlation 
is present, with more massive objects having on average the 
larger cross sections. The scatter on the correlation is however 
considerably large, being up to one order of magnitude both in 
mass and cross section. This scatter is due to the combination 
of many different effects, including cluster ellipticity, substruc- 
tures, and alignment with the line of sight. We investigate more 
on this scatter further below. 

In Figure |2] we also reported the expected trend with mass 
and redshift of the cross section for giant arcs produced by a 
Singular Isothermal Sphere (SIS henceforth) lens model act- 
ingon circular, point-like sources (see Naravan & Bartelmann 
19991 for a thorough description of this simple lens model). For 
arcs with length-to-width ratio larger than Rq, this cross section 
can be written as 



o-Ro = 



327r'Dl{zL,Z5) 



Rl + l 



(Rl - 1)^ 



(6) 



where ctdm is the velocity dispersion of the SIS density profile, 
that we identify with the average velocity dispersion of dark 



matter particles within the cluster. We adopted the scaling rela- 
tion of Evrard et alj (120081) in order to relate ctdm to the mass 



of clusters, obtaining that 



O-Ro = 



32n^Dl{zL,Zs) /?g + 1 
c4 (Ri - 1)2 



CTn 



h(ZL)M 



lO'^Mo 



4a 



(7) 



In the previous equation we used a - 0.336 and cr^Mfi - 
1.083 X 10"* km s"', while /i(zl) = H(zO/Ho. The amplitude 
of the cross sections for the SIS models is almost always larger 
than the best fit we have found, and this is particularly evident 
at low redshift. Had we adopted elliptical extended sources for 
the SIS lenses as well, the relative cross sections would have 
been increased further, thus worsening the disagreement. In 
part this can be explained by the fact that what is reported in 
Figure |2] is the average of the lensing cross sections over three 
orthogonal projections for the MareNostrum Universe clus- 
ters. Since it is very unlikely that all the three cross sections 
for each cluster are different from zero, this would downweight 
the average cross section with respect to the case in which the 
three cross sections have similar magnitudes, as is the case for 
SIS lenses (in which obviously ctk^ is the same for all projec- 
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Fig. 4. The correlation between the logarithm of the cross section for giant arcs (y - log(a"7 5 li^/Mpc^)) and the logarithm of the 
bolometric X-ray luminosity (in units of lO"*"* erg s"', ^ = log[Lx s/(10'*'*erg)]) for various redshift snapshots of the MareNostrum 
Universe simulation. The red lines show the best linear fits, while the green dashed lines show the self-similar expectation detiled 
in the text. 



tions). Even this however cannot remove the whole disagree- 
ment, especially at low redshift and especially if we consider 
source ellipticity and finite size. Another and more physical 
effect is that, due to the unrealistically steep density profile, 
the SIS cross section at low redshifts is expected to be always 
larger than the cross section of numerically simulated clus- 
ters , that tend instead to follow more closely an internally flat- 
ter (INavarro et al.l 19961 N FW henceforth) density pr ofile. This 
fact was already noted by iMeneghetti et al.l (l2003bl) . where it 
was also shown that at high redshifts the steepness of SIS lenses 
compensates quite well the inadequacy of spherical density 
profiles to reproduce the lensing properties of realistic halos, 
giving rise to cross sections that are in good agreement with 
those of the latter. This fact also is confirmed in our Figure |2] 

Despite the different normalizations, the slopes of the best 
fit X - y relations and those derived for the SIS models are in- 
triguingly similar (despite the latter being slightly flatter), es- 
pecially in light of the large scatter displayed in Figure |2l This 
impUes the interesting fact that, notwithstanding the spread that 
is introduced in the cluster population by differences in individ- 
ual formation histories (mergers, substructures, ellipticity, etc.). 



the simple "self-similar" scaling o-Rg cc R^ oc cr^j^ oc M'^^^ 
(where Re is the Einstein radius of the lens) can be considered 
a fair representation of average strong lensing clusters. 

As noted above the two parameters of the linear fit change 
with the lens redshift, however we did not find any discernible 
trend, with both the a{zL) and b(zL) functions being rather flat 
and oscillating around almost fixed values. An exception to this 
is given by the very low or very high-redshift snapshots, where 
the fluctuations become very large. However in those snapshots 
very few cluster with non-vanishing cross sections are avail- 
able, hence the resulting fit is likely not to be trusted. We come 
back to this issue again in Section l4~2l 

Since the scatter around the best fit relations in Figure |2]is 
substantial, it is necessary to quantify it. Hence, for each strong 
lensing cluster present in each redshift snapshot we computed 
the quantity y - y^, where y^ is the logarithm of the cross sec- 
tion given by the best fit relation at that redshift and for a fixed 
logarithmic mass, while y is the logarithm of the actual cross 
section. In Figure|3]we show how, for each of the redshift snap- 
shots shown in Figure |2] these quantities are distributed. We 
attempted at first to fit these distributions with Gaussians, with 
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Fig. 5. The distribution of the logarithm of cross sections for giant arcs (y = log(a-7 5 /i^/Mpc^)) around the best fit relation at a 
given X-ray luminosity (^ = log[Lx s/(10'*'*erg)]). The green dashed line represents the best Gaussian fit, while the solid red lines 
are the best non-Gaussian corrections. 



the resulting best fitting curves shown by green dashed curves 
in Figure [3] The fit is fair but not exceptional, given the rather 
irregular aspect of the distributions, especially for the high and 
low-redshift snapshots. Therefore, we added one free param- 
eter representing a suitably small deviation from a Gaussian 
distribution. In order for this deviation to be related to the 
moments of the distribution, we adop ted the functional form 
suggested by the Edgeworth expansion JBlinnikov & Moessnei 
1998|;ILo Verde et al.ll2008l) . namely 



p(y - yx) 



1 



^/27TC 



■ exp 



(y-y.f 



2cr2 



1 +o-Hi,(y-yo)- 



(8) 



where H^iu) is the Hermite polynomial of third degree, that 
reads 



Hj,{u) — u — 3m. 



(9) 



The parameter ^3 = Sj,{zC) is the normalized skewness of the 
distribution, and it has been left as the second free parameter 
Obviously if 5 3 = we fold back to the Gaussian case where 
the only free parameter is cr = cr{zi) (not to be confused with 
the lensing cross section). 



In Figure [3] the red solid curves show the results of this 
kind of non-Gaussian fit. As can be seen, Eq. ([8]) provides a 
slight improvement over the simply Gaussian fit, especially 
for low and high-redshift snapshots. It is likely that the fit 
could be improved further by introducing a third free param- 
eter (that in the Edgeworth expansion would represent the nor- 
malized kurtosis). Also, there are several snapshots (e.g., those 
at z = 0.56 and z - 0.66) that are showing traces of bimodal- 
ity, for which a fit with the superposition of two curves would 
be preferable. As we discuss in Section |4|2] the two parameter 
fits that are shown in Figure [3] are not precise enough to accu- 
rately reproduce the strong lensing optical depth measured in 
the MareNostrum Simulation. Since however we have found 
that the inclusion of the scatter is not necessary for this scope, 
we prefer to add no further complication and stick to the two 
parameter fits. Besides, while the redshift evolution of crizi) is 
rather flat around an average of ~ 0.3, the redshift evolution 
of the parameter 5 3 (zl) is much more noisy than for fl(zL) and 
^(zl), showing large fluctuations between neighboring redshift 
snapshots. This fact suggests that with the introduction of the 
skewness as an additional fit parameter we are in part overfit- 
ting the noise of measured distributions. The addition of further 
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parameters would worsen this situation even more, thus being 
of little practical use. 

4.1.2. Luminosity 

Since the MareNostrum Universe simulation includes adia- 
batic gas, we could consider the coiTelations of the virial mass 
and the strong lensing cross section with the bolometric X-ray 
luminosity, the latter being a quantity that is more directly re- 
lated to observable features. To that end, we introduced a new 
variable ^ = log Lx s/ (lO"^'*ergj , where Lx is the bolometric 
X-ray luminosity. In Figure |4] we report the relation between y 
and ^, while the ^ - x relation has been already discussed in 
Paper I. 

In the same Figure we also report the best linear fit (red 
solid lines) and the y - ^ relation that can be deduced by com- 
bining the SIS scaling derived above for the y- x relation with 
self-similar arguments (green dashed lines), according to the 
following line of reasoning. Let us assume that dark matter ha- 
les in the MareNostrum Universe simulation follow on av- 
erage a NFW density profile. Self-similar scaling arguments 
imply that the bolometric luminosity of a plasma in hydrostatic 
equilibrium in the potential well of a dark matter halo with s uch 
a density profile follows JEke et al.lll998HFedeli et aTllioOTl) 
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where zl is the redshift of the lensing cluster at hand, M is its 
virial mass, and c is its concentration. The function F(c) reads 



F(c) = ln(l-Hc)- 



1 +c 



(11) 



Since even for the most massive clusters the concentration is 
at least of a few, we have that F(c) ^ ln(c) - 1, and this de- 
pendence, being very weak, can be neglected. We further recall 
that, roughly, c oc M""' thus, neglecting the redshift depen- 
dence of the concentration, overall we obtain the approximate 
scaling Lx(M,zl) °^ Mh{z]y^^- In combination with Eq. ^ 
this finally implies 
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In order to fix the normalization we simply required that the 
bolometric X-ray luminosity and cross sections of the most 
luminous cluster in the simulation are reproduced. We would 
like to stress that the procedure we followed for deriving of the 
above scaling relations is likely oversimplified. Here we did not 
want to perform a completely rigorous derivation, just to make 
the case for the fact that simple scaling arguments can produce 
coiTelations that are rather similar in slope to those obsei-ved in 
the MareNostrum Universe. 

Because of the increase in average X-ray luminosity with 
increasing mass and because of the trend reported in Figure |2] 
it is expected that the lensing cross section also increases with 
the X-ray emission. This is indeed observed in Figure |4] and 
the correlation is visibly tighter than the y — x correlation. In 



order to better appreciate this, in Figure |5] we report the distri- 
bution of logarithmic cross sections around the best fit relation 
for a fixed X-ray luminosity. As a matter of fact, the logarith- 
mic slope of the y - ^ relation is relatively similar to that of 
the y-x relation (in agreement with the self-similar scaling de- 
tailed above), but the scatter around the latter is systematically 
~ 25% smaller than around the former. Also, as we already 
found for the mass-cross section relation, the introduction of 
a non-Gaussian part to the functional form used for fitting the 
distribution around the best suitable luminosity-cross section 
relation does not substantially alter the fit with respect to the 
perfectly Gaussian form. 

The fact that the scatter around the best fit relation is de- 
creased suggests that there should be a more fundamental re- 
lation of the cross section for giant arcs with the bolometric 
X-ray luminosity, rather than with the virial mass. This can be 
somewhat expected from several qualitative considerations: a) 
relaxed clusters are expected to have large gas and dark matter 
central densities. This would increase both the X-ray emissiv- 
ity, that is proportional to the squ are of the gas den sity, and 
the extension of the critical curves (iFedeh et al.ll2007l) : b) clus- 
ter mergers are expected to compress and heat up the intra- 
cluster medium, hence increasing the X-ray luminosity. At the 
same time, interactions s hould increase th e convergence and 
shear fields of the cluster (ITorri et al.ll2004l) . thus enhancing its 
lensing efliciency; c) in general, the X-ray luminosity is par- 
ticularly sensitive to the physical processes taking place in the 
inner region of galaxy clusters, and so is the strong lensing effi- 
ciency. The virial mass on the other hand is much less sensitive 
to that. This also implies that the correlation between the lens- 
ing cross section and the mass contained within smaller radii 
(e.g., M5()o,M25oo, etc.), should be tighter than with the virial 
mass. 

A closer inspection of Figure |4] shows that virtually all the 
clusters expected to have at least one giant arc (0-75 > 10"^ 
Mpc^/!"^) have a bolometric X-ray luminosity that is larger 
than ~ 6 X lO"*^ erg s"'. The converse is in general not true, 
in the sense that many clusters exist with a high luminosity that 
make rather poor strong lenses. Hence, strong lensing clusters 
would be good tracers of centrally dense intra-cluster medium, 
while the converse would not hold in general. 

Finally, we tried to recompute the scaling relations between 
X and ^ and the logarithm of the cross section y by taking 
into account only relaxed clusters. Relaxation was quantified 
through the parameter yS introduced in Paper I, that expresses 
deviations from virial equilibrium: negative p implies systems 
dominated by the kinetic energy, positive yS by the potential en- 
ergy, while 13 -Q implies perfect virial equilibrium. We recom- 
puted the above scaling relations by using only clusters with 
|y6| < 0.3, and found no reduction in the scatter about the best 
fit X - y and ^ - y relations, nor a significant shift of the quar- 
tiles of the distributions. What happen instead is just a removal 
of the highest mass/luminosity objects, which is expected since 
those are the least relaxed on average. This is an indication that, 
intriguingly, the bulk of this scatter should not be originated by 
the dynamical activity of clusters, rather by triaxiality and con- 
centration distributions. The latter, particularly, is expected to 
have a major impact since the cross section is remarkably sen- 
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sitive to it, and the concentration itself has a substantial intrin- 
sic scatter. Because of the relatively low resolution of clusters 
in the simulation, we could compute only concentrations for 
staked cluster samples, not for individual objects. Hence we 
could not asses this issue in detail. 

4.2. Optical depths 

4.2.1 . Synthetic optical deptii 

In Figure|6]we show the differential optical depth as a function 
of lens redshift. The differential optical depth is the integrand of 
Eq. ^, and represents the contribution to the total optical depth 
that is given by structures at different lens redshifts. As a conse- 
quence, the integral under the black solid curve in Figure|6]rep- 
resents the total optical depth, that for sources at zs =2 amount 
to T7.5(zs = 2) ^ 1.59 X 10 ''. The differential optical depth has 
the generic shape that is expected, namely dropping at low and 
high redshift due to the geometric suppression of the lensing 
efficiency. The total optical depth is subs tantially lower tha n 
that reported for a erg = 0.9 cosmology in lFedeli et al.l (12008'), 
which is to be expected since there the authors used a complete 
source redshift distribution, while here we adopted sources at 
a fixed redshift. By multiplying the total optical depth with the 
number counts of source galaxies at zs -2 it is possible to infer 
the abundance of arcs produced by those sources that are visi- 
ble below some given magnitude threshold (see the discussion 
is Section[2]i. 

The existence of a relation between the virial mass of 
galaxy clusters and their cross sections for giant arcs (Section 
14. Il l implies that we could in principle construct a synthetic re- 
alization of the optical depth in the following way. For each 
step in redshift we randomly extract a large number of masses, 
associate a cross section to each mass through the jc - y re- 
lation described above, suitably randomized according to the 
scatter around the best fit relation, and use these cross sec- 
tions for computing the differential optical depth at the fixed 
redshift. This kind of procedure would allow one to compute 
strong lensing optical depths based only on the scaling relations 
that we have found and without expensive calculations involv- 
ing ray-tracing simulations and cross section evaluations. 

In Figure [T] we show the redshift evolution of the param- 
eters a(zL) and ^(zl) defined in the previous Section 14.11 As 
announced, both of them oscillate around a roughly constant 
value. We plot in the same Figure also the average values, com- 
puted over the redshift range [0.1,1.1], since outside this range 
very little lensing clusters are present per redshift snapshot, and 
the result of the fitting procedure is less reliable. These two av- 
erage values amount to a = 1.529 and b = -26.08. Although 
the procedure described above would make use of the specific 
vales a(zL) and b{zL) at each redshift zl, this would not be of 
very practical use. Therefore we decided to adopt the same val- 
ues of the two parameters a and b at each redshift snapshot, 
equal to the respective average values. At first, we tried to ig- 
nore the scatter around the best fit y - x relation, i.e., we as- 
signed at each mass the exact value that is prescribed by the 
linear logarithmic fit with the redshift averaged values. Another 
factor that needs to be taken into account is the smallest mass 
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Fig. 7. The redshift evolution of the two parameters defining the 
best fit x-y relation (see the text for details), with the respective 
standard errors. The black and red solid horizontal lines repre- 
sent the respective redshift averages, while the green horizontal 
line represents the slope expected for a SIS lens model. 

MminizL) that, at each redshift step, enters in the computation 
of the differential optical depth (Mi in Eq. |2]i. As it turns out 
this point is of fundamental importance, since the optical depth 
tends to be dominated by the smallest clusters capable of pro- 
ducing a non-vanishing cross section, being by far the more 
abundant. Instead of adopting as such lower limit the small- 
est mass that at each redshift produces a non-vanishing cross 
section in the MareNostrum Universe we let the redshift evo- 
lution of M,nin(zL) be 



s(Zh.p) 



s{zl) 



(13) 



We adopted this procedure for a twofold reason. Firstly, choos- 
ing a different minimum mass at each redshift snapshot is not 
practical, exactly as it was not practical to use a different pair 
of best fitting parameters at each redshift. Second, it is likely 
that this minimum mass is siumlation-dependent, or in different 
words, that rerunning the same simulation several times would 
produce a distribution of minimum masses at a given redshift, 
due to the randomness of the cluster formation histories. 

In Eq. ( fTlT l the function eizh) is the lensing efficiency that 
has been introduced in Section [3] zl.p is some pivotal lens red- 
shift and Mmin.p is the corresponding minimum lensing mass. 
We adopted as pivotal the redshift at which the differential op- 
tical depth is maximal, zl.p - 0.53, that also coincides with the 
broad maximum in the lensing efficiency reported in Figure [1] 
At this redshift the smallest mass with non- vanishing cross sec- 
tion is M, 



10 Mq/; . Ideally, selecting a suitable value 
of q for each lens redshift would allow to match the synthetic 
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Fig. 6. Left panel. The differential optical depth for giant arcs obtained by all clusters in the MareNostrum Universe as a 
function of lens redshift (black solid line). The red dashed line represents the synthetic differential optical depth obtained with 
the procedure outlined in the text {q - 3.5), while the shaded area represents the effect of a 30% uncertainty on the minimum 
mass at the pivotal redshift (see the text for details). Right panel. Same as the left panel but with the blue and green dashed lines 
now representing the synthetic optical depths computed with values of q different from 3.5, as labeled. 



optical depth, since one would encapsulate in this free param- 
eter all the remaining uncertainty that has not been taken into 
account, i.e., the deviation of the parameters ^(zl) and ^(zl) 
from their redshift averages, the scatter around the best fit x - y 
relation, etc. In practice, we can choose a single value of q that 
gives an overall good fit to the true optical depth. 

The "synthetic" optical depth resulting from this procedure 
is represented in Figure |6] together with the effect of allow- 
ing some uncertainty on the value of Mmin,p and changing q 
within some reasonable interval. The deviation of the synthetic 
optical depth from the true one for, e.g., q - 3.5 might seem 
quite significant. Particularly, the true optical depth appears to 
be quite overestimated at zl = 0.2 - 0.3 due to the fact that 
both the parameters a{zi) and bizi) are below the average at 
those redshifts (see Figure |7]). However it should be recalled 
that the Figure has logarithmic scale, and as a matter of fact the 
synthetic optical depth with ^ = 3.5 differs from the true one 
only by about ~ 10 - 15%, which is a truly negligible amount 
given the large uncertainties involved in the estimation of arc 
abundances. Increasing the value of q to, e.g., 4.5 reduces the 
discrepancy on the total optical depth to only ~ 5%, however in 
this case the high-redshift tail of the differential optical depth is 
rather underestimated. Based on these numbers and on Figure 
|6]we suggest to adopt a value of ^ ^ 3.5 as a fair compromise. 

Also shown in Figure |6] is the effect of uncertainty in the 
value of Mmin.p. As a matter of fact, a substantial fluctuation 
on the minimum lensing mass at redshifts close to the pivotal 
one is found, due to the fact that in such a boundary situation 
small variations in the internal structure of the lens are enough 
to scatter the lens itself above or below the threshold for cross 



section production. In Figure |6] we show as an example the ef- 
fect of a 30% uncertainty on Mn,in,p, which can easily shift the 
total optical depth of ~ 50%. The effect of changing the value 
of q within a reasonable interval is relatively similar in ampli- 
tude to the effect of Mmin.p fluctuations. 

We have to observe that the reasonable way in which the 
synthetic optical depth approximates the true one is heavily de- 
pendent on the values of a and b that we assume. Particularly, 
to produce variations to the synthetic optical depth with mag- 
nitude similar to that reported in the left panel of Figure |6] it is 
sufficient to modify the normalization or the slope of the mass- 
cross section relation at the percent level. This is way smaller 
than the fluctuations and errors that are reported in Figure I?] 
For instance, adopting the slope of the mass-cross section re- 
lation that is found for SIS lens model, instead of the redshift 
average that we adopted, would produce a gross underestimate 
of the true optical depth. This argument can be turned around 
by saying that the synthetic optical depth is a fair representa- 
tion of the true one only if the normalization and the slope of 
the mass-cross section relation adopted stay very close to their 
redshift-average values computed here. 

Lastly, we would like to comment on the fact that in the 
above procedure for computing synthetic optical depth we ne- 
glected the role of the scatter around the best fit y - x relation. 
As we mentioned in Section 14.11 the fits that we performed 
to the distributions of y - y^ both with a Gaussian and a non- 
Gaussian functions are not precise enough for the purpose of 
increasing the precision of the synthetic optical depth. This can 
be understood as follows: as it appears evident from a close in- 
spection of Figure[3l neither the Gaussian nor the non-Gaussian 
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Fig. 8. The difference between the true distribution around the best fit mass-cross section relation and the best Gaussian (blue 
shaded histograms with black contours) and non-Gaussian (magenta shaded histograms) fits. Results for several different redshift 
snapshots are shown, as labeled in the plot. 



functional forms fully appreciate the asymmetry in the scatter 
of cross sections around the best fit. Namely, the distributions 
of y - y.-c tend to be substantially skewed toward negative val- 
ues as compared to the best fits. As a consequence, the intro- 
duction of the scatter around the best fit actually worsen the 
agreement of the synthetic optical depth to the true one. This 
fact can be better appreciated by analyzing Figure[8] There, the 
scatter around the best fit mass-cross section relation is shown 
once the best Gaussian and non-Gaussian fits were removed. It 
is obvious that the best Gaussian fit is substantially more posi- 
tively skewed with respect to the true distribution, especially at 
high redshift. The same remains true for the non-Gaussian fit, 
although to a lesser extent. This means that, including the scat- 
ter in the computation of the synthetic optical depth effectively 
equals at assigning higher cross sections to a given mass with 
respect to reality. As a result the synthetic optical depth tends 
to be a substantial overestimate of the true one if the scatter is 
included 

In order to solve this, a better fit to the distributions of y-y^ 
would be required, maybe with the superposition of two curves. 
This however would increase the number of free parameters 
from two to four at least. We chose not to pursue this route, 



since it would add substantial complication to a model that al- 
ready works quite acceptably. As a matter of fact, given the 
very simplistic nature of the underlying assumptions (one best 
fit relation valid for all redshifts, absence of scatter, etc.) it is 
very remarkable that the true optical depth can be reproduced at 
this level with just one adjustable parameter, q. This very sim- 
ple model for computing synthetic optical depths can be very 
helpful for, e.g., evaluate the contribution to arc abundance that 
is given by structures in various redshift ranges without per- 
forming any actual strong lensing calculation. Summarizing, 
in evaluating synthetic optical depths for sources at zs =2 we 
suggest to stick to the precise values of a and h that we adopted, 
and to use a value of q between 3.5 and 4.5, with the threshold 
mass at the pivot redshift equaling IO'^^Mq/z"'. 

4.2.2. The role of dynamical activity 

In Figure |9] we show the differential optical depth of the 
MareNostrum Universe obtained by excising those clusters 
that are not relaxed enough. The relaxation of a cluster was es- 
tablished vi a the viri al equilibrium parameter, defined in Paper 
1 (see also IShaw et al. 2006) . As can be seen the exclusion 
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Fig. 9. The differential optical depth in the MareNostrum 
Universe computed for all clusters in the simulation (black 
solid line) and by including only those clusters whose (3 param- 
eter is included in some neighborhood of zero (dashed colored 
Unes), as labeled in the plot. 



of clusters with a value of p significantly different from zero 
causes an erosion of the high-redshift tail of the differential op- 
tical depth. Particularly, when only clusters with \[i\ < 0.6 are 
included, the differential optical depth peaks at Zh ~ 0.4 instead 
of zl ~ 0.55, and the total optical depth is reduced by a factor 
of 2. In other words, half of the number of giant arcs produced 
by sources at zs =2 are expected to be found inside clusters 
that are substantially unrelaxed. A significant contribution to 
arc statistics from clusters at zl > 0.8 can be expected only by 
including objects with ^| < 1. 

This trend of the optical depth with the virial equilibrium 
of clusters can be easily understood by looking at the redshift 
evolution of/3 that has been presented in Paper 1. Namely, while 
at z ~ the average /3 is quite close to zero, at high redshifts it 
tends to become strongly negative. This means that, removing 
clusters with \/3\ very different from zero is equivalent at remov- 
ing high-redshift clusters, and hence the drop in the differential 
optical depth at large zl follows. Intriguingly, the redshift dis- 
tribution of strong lensing clusters carries information on the 
dynamical state of the deflectors: selecting clusters that are ca- 
pable of producing giant gravitational arcs at zl ^ 0.5 we are 
automatically selecting objects that are relatively close to the 
virial equilibrium. Vice versa, lensing clusters at high redshift 
are very likely to be dynamically active. 

We wish to interpret this result in light of the incidence 
of gravitationa l arcs in high-red shift clusters found by, e.g ., 
Gladders et al] (l2003b (see also [zTritskv & Gonzalez! I2OO3I) . 
Particularly, Gladders et al. (2003) analyzed the Red-sequence 
Cluster Survey (RCS, Gladders & Yee 2001) . finding eight 



clusters at zl ^ 0.64 with prominent gravitational arcs, for a to- 
tal of 11 strong lensing features, including tentative detections. 
By using Figure |9] we can estimate that ~ 77% of the optical 
depth at zl ^ 0.64 is produced by clusters with \/3\ > 0.6, which 
effectively means /3 < -0.6 since the values of /3 are promi- 
nently negative, especially at high redshift (see Paper I). Thus, 
it is ex pected that ~ 8 - 9 of the 1 1 giant arcs of lGladders et al 



( l2003i) are pro duced by substantially unrelaxed clusters. As a 
matter of fact, iGladders et al.l (l2003h themselves argued that 
the high-redshift lensing clusters they found are somewhat en- 
hanced in their lensing efficiency by some physical process, 
possibly the presence of substructures and dynamical activity. 

Browsing the literature in order to confirm the fraction of 
strong lensing clusters that are dynamically unrelaxed gives 
inconclusive results. For instance, RCS 2319.9-H0038 is part 
of a supercluster and is likely to b e unrelaxed, as well as 
RCS 2156 .7-0448 JHicks et alj |2"008l). RCS 1620.2-^2929 was 
shown by Gilbank et al.l (l2007h to have an excess velocity dis- 
persion c ompared to what wou ld be expected by its optical 
richness (lYee & Ellingsonll2003h . suggesting signs of dynam- 
ical activity. O n the other hand RCS 0224.5-0002 is Ukely a 
relaxed object JHicks et al.l2007l) . totaling to 4 giant arcs found 
in dynamically active clusters versus 2 found in relaxed clus- 
ters. All the other RCS clusters used bv lGladders et al.l (l2003b 
do not have information about the dynamical state. Obviously 
we are not able to draw any conclusion from this, but just want 
to stress that further study of the dynamical activity of strong 
lensing clusters, especially at high redshift, would be valuable 
in order to better characterize the cluster lens population. 

4.2.3. X-ray selection effects 

It is finally interesting to investigate how the optical depth 
changes if clusters are selected with X-ray luminosity. 
This is particularly worthy since past observational arc 
statistics s tudies focuse d mainly on X-ray selected clusters 
(.Gioia & Luppinol 1 1994c iLe Fevre et alJ 1 1994c iLuppino et al 



19991) . In Figure [To] we show the fraction of the total optical 



depth for giant arcs that is contributed by clusters with lumi- 
nosity above a given threshold. Under the assumption that all 
sources lie at zs =2 this corresponds to the fraction of gravi- 
tational arcs that are found inside luminous clusters. Two lines 
are shown in Figure [TOl one referring to the bolometric X-ray 
luminosity (labeled Lx), the other to the luminosity computed 
in the [0.5,2] keV energy band (labeled Lb). The latter was 
estimated from t he integrated luminos ity by using the bolomet- 
ric correction of iBorgani et al.l (119991) . For that we needed the 
global cluster temperature, which we estimated from the red- 
shift and the virial mass of the object by using the virial rela- 
tion 
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where Av(zl) is the virial overdensity of collapsed top-hat den- 
sity perturbations at redshift zl and the parameter y defines the 
ratio of the kinetic energy per unit mass of the dark matter to 
that of the gas. For simplicity we adopted y - I henceforth 
(lEvrard et al. 2008) . 
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Since no clusters with Lx ^ 10 erg s have non- 
vanishing cross section (see Figure|4]l, the totality of the strong 
lensing optical depth is produced by clusters with bolometric 
luminosity above that value, as can be seen in Figure [10] In 
other words, if one considers a complete cluster sample con- 
taining all the objects with Lx > 10'*'* erg s ', then this sample 
would include all the strong lensing clusters in the surveyed 
area. While several clusters with Lx > 10"*^ erg s"' exist, their 
contribution to the total optical depth is small, since they are 
likely to be very massive and hence extremely rare. As a mat- 
ter of fact. Figure [To] shows that such objects contribute only 
~ 20% of the total giant arc abundance. If one considers the 
band luminosity, clusters with Lb > 10^"* erg s"' also contribute 
more than 90% of the total arc abundance. However, the drop 
of the optical depth with luminosity is somewhat steeper than 
before, with clusters having Lb > 3 x 10'*'* erg s"' contributing 
only ~ 20% to the total optical depth. 

These results are particularly relevant in order to estab- 
lish the reliability of observational studies directed toward 
arc statistics. As a classic example (and witho ut the pretense 



of bein g rigorous), we considered the work of iLe Fevre et al. 



(119941) . where the authors performed a search for giant arcs 
inside clusters extracted from the Einst ein Medium Sensitivity 
Survey (EMSS. lGioia & Luppinol 19941) by having luminosities 



in the [0.3, 3.5] keV energy band larger than 4 x 10'*'* erg s '. 
Although we did not compute the luminosities in that energy 
band for all of our simulated clusters , the relative result will lie 
somewhere in between the red and black lines in Figure[TO] im- 
plying that at least ~ 20% of the arcs were missed in that sam- 
ple. We estimated the true luminosity in the [0.3, 3.5] keV en- 
ergy band only for a subsample of the MareNostrum Universe 
lensing clusters, and found that this fraction should be about 
30 - 40%. Hence, these observational estimates of arc abun- 
dances would have to be corrected for this factor, thus wors- 



ting the disagreement with theoretical estimates (IFedeli et al 



12008 *). Nevertheless, such a correction is likely to be small 
compared to other uncertainties concerning theoretical and ob- 
servational estimates of arc statistics. 

This kind of stud y is somewh a t simi lar to the analysis 
performed in .Fedeli & BartelmannI (l2007l) . where the semi- 
analytic optical depth actually observed in flux limited X-ray 
cluster samples was detailed. However, the two works cannot 
be directly compared because there a full account for instru- 
mental effects was pre sented that is missing here and a dififerent 
plasma model is used JRavmond & Smith]] 1977b . 



5. Summary and discussion 

In this work we investigated the optical depth emerging from 
the strong lensing properties of a very large set of numeri- 
cal galaxy clusters extracted from the MareNostrum Universe 
cosmological simulation. The simulation has the same box- 
size of the Millennium run fepringell ]2005h . and includes 
an equal number of particles for both dark matter and adia- 
batic gas. The cosmological parameters used in the simula- 
tion are in agreement with the WMAP-1 year data release, 
and in particular the normalisation of the matter power spec- 
trum is cr^ = 0.9. Although this is higher than the results 
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Fig. 10. The fraction of the total optical depth that is con- 
tributed by clusters of the MareNostrum Universe with log- 
arithm of the bolometric X-ray luminosity (in units of 10'*'* erg 
s"', ^ = log[Lx s/(10'*'*erg)]) larger than the value reported on 
the abscissa (filled circles connected by the black line). The 
filled triangles connected by the red line show the same quan- 
tity evaluated for the luminosity in the energy band [0.5, 2] keV. 



of t he WMAP-5 year dat a analysis dKomatsu et al.]]2009l see 



also Komatsu et al 



2010), the two are still consistent at 95% 
Confidence Level. Additionally, studies of the present abun- 
dance of X-ray clusters indicate a somew hat higher normal- 
ization than found by the W MAP team ( Yepes et al.l 12007 : 
Wen et al.ll2010HWoitak & Lokas 2010) . 

Either way, a lower normaUzation of the matter power spec- 
trum would imply a later formation of cosmic structure, with 
consequent lower abundance of massive galaxy clusters and 
a lower average concentration there of. Both of these effects 
would act to reduce the optical depth (F edeli et al . 2008), how- 
ever, we expect the relative contributions to the arc abundance 
given by unrelaxed/luminous clusters to be mildly aff'ected. 
Additionally, a lower dark matter halo concentration might be 
compensated by the introduction of gas cooling, that would 
increase the amount of mass flowing to the central region of 
galaxy clusters. 

Before summing up, there is one point of our analysis that 
we would Uke to stress. The population of simulated dark mat- 
ter halos capable of producing giant arcs has a lower mass limit 
due to the fact that when caustic curves become too small com- 
pared to the average source size, the related images can only 
very rarely be strongly distorted. For this reason isolated mas- 
sive galaxies, that yet can produce arcs in particularly favor- 
able circumstances, did not enter in the computation of the op- 
tical depths. Similarly, the simulation cannot realistically ac- 
count for the formation of individual galaxies inside clusters 
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due to the lack of baryon cooling, and therefore giant arcs that 
are mainly produced by those galaxies are also not included 
in our computations. However we argue for the influence of 
these issues on the optical depths to be small, since the ef- 
fect of individual galaxies on the stron g lensing cross sections 
of clusters has already b een estimated dMeneghetti et alJ l2Q00; 
Wambsganss et al.ll2008 ) to be mild, and the probability of iso- 
lated galaxies for producing giant arcs is much too small for 
it being compensated by their larger abundance. Besides, real 
strong cluster lensing observations do not focus, by definition, 
on lensing by isolated galaxies. Nevertheless, the aforemen- 
tioned limitations of our analysis are to be kept in mind. 

The main conclusions that we reached in this paper can be 
summarized as follows. 

• We quantified the correlation between the virial mass of 
galaxy clusters and the cross section for giant arcs. This 
relation is approximately linear in logarithm with a slope 
around ~ 1 .5 and accounts for the expected fact that more 
massive objects have an higher probability of producing 
strong lensing events. The parameters of the best fit are 
roughly constant with redshift, although significant oscil- 
lations around the mean are displayed. 

• A similar correlation (with a similar slope) was found be- 
tween the bolometric X-ray luminosity and the lensing effi- 
ciency. This correlation appears to the tighter than the pre- 
vious one, reflecting the fact that both X-ray emissivity and 
lensing cross sections are particularly sensitive to the pro- 
cesses occurring at the very center of galaxy clusters. 

• Both the correlations described above display a substantial 
scatter. In particular, the scatter around the best fit mass- 
cross section relation can be up to one order of magnitude. 
We found that the distributions around the best fit logarith- 
mic relations are approximately Gaussian, although signif- 
icant deviations from Gaussianity might be found. 

• Using the best fit correlations between mass and lensing ef- 
ficiency we managed to construct a synthetic optical depth 
that with only one adjustable parameter manages to repro- 
duce acceptably well the true differential optical depth. 

• We computed the contribution to the total optical depth that 
is given by relaxed and unrelaxed structures, where relax- 
ation was defined as deviation from virial equilibrium. We 
found that unrelaxed structures mainly contribute to the 
high-redshift part of the optical depth, in agreement with 
the standard paradigm for structure formation. 

• Similarly, we evaluated the contribution to the total arc 
abundance given by structures of various X-ray luminos- 
ity, thus quantifying the bias present in arc statistic studies 
based on X-ray selected clusters. 



Recently iHoresh et al.l (120101) have shown that X-ray se- 
lected clusters display a substantially higher giant arc incidence 
than optically selected clusters with comparable optical lumi- 
nosity. They conclude that X-ray observables trace much larger 
mass concentrations than optical ones. Our finding of a rela- 
tively tight correlation between the lensing efficiency and the 
bolometric X-ray luminosity go es in the same direc tion, and 
can help interpret the findings of lHoresh et alJ (l2010l) . 



We believe the large set of strong lensing cross sections 
computed in this work to have a wide range of future applica- 
tions beyond those exploited here and in Paper 1, since it covers 
an abundant sample of realistically simulated galaxy clusters. 
We are currently planning to use cluster observables like the X- 
ray luminosity to construct realistic past light cones, in order to 
better evaluate possible selection effects and to better compare 
the lens population to the overall cluster population. 
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